Superconducting transport through a vibrating molecule 
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Nonequilibrium electronic transport through a molecular level weakly coupled to a single coherent 
phonon/ vibration mode has been studied for superconducting leads. The Keldysh Green function 
formalism is used to compute the current for the entire bias voltage range. In the subgap regime, 
Multiple Andreev Reflection (MAR) processes accompanied by phonon emission cause rich structure 
near the onset of MAR channels, including an even-odd parity effect that can be interpreted in 
terms of an inelastic MAR ladder picture. Thereby we establish a connection between the Keldysh 
formalism and the Landauer scattering approach for inelastic MAR. 

PACS numbers: 74.78.Na, 74.45.+C, 74.50.+r 



I. INTRODUCTION 

One of the primary goals in the field of molecular elec- 
tronics is to understand quantum transport through indi- 
vidual nanoscale objects, such as molecules, short carbon 
nanotubes, or DNA; for reviews, see Refs. 0,111!. An 
important difference to conventional mesoscopic trans- 
port through quantum dots or granular islands arises be- 
cause molecules can have intrinsic vibrational degrees of 
freedom ('phonons') that may give rise to Franck-Condon 
sidebands or phonon blockade in electronic transport. 
Such features have been studied theoretically in many 
recent papers 0, El E3, 0, El- Molecular electronics is 
particularly exciting because it is in principle possible 
to contact molecules by leads of different nature. Here 
we discuss how nonequilibrium transport is affected by a 
coherent phonon mode coupled to the molecular charge 
for the case of (s-wave BCS) superconducting leads. 
Molecules connected to superconductors promise a rich 
terrain of exploration that allows for new spectroscopic 
tools (probing molecular properties), potentially useful 
applications, and interesting fundamental physics. First 
ex peri mental results have appeared for carbon nanotubes 
013111 Elll and metallofullerenes Q. 

So far transport through molecules has been theoret- 
ically studied only for normal leads, either using rate 
equations or (in the quantum-coherent regime Vperturba- 
tion theory in the electron-phonon coupling A • For su- 
perconducting leads and large transmission through the 
molecule, subgap transport is ruled by MAR processes 
[Tol ITt| . These have been extensively studied for point 
contacts Il7l Il8l Il9l and for junctions containing a reso- 
nant level |20l l2lj . In this paper, we provide a theoret- 
ical framework to include vibrations into superconduct- 
ing transport through a resonant molecular level. We 
focus on the most interesting quantum-coherent low tem- 
perature (T) limit with high (bare) transmission, where 
Coulomb charging effects are largely wiped out. There- 
fore the Coulomb interaction on the molecule will be ne- 
glected here. We compute the d.c. current basically for 



the full bias (V) range within a Keldysh Green function 
scheme valid for small A but arbitrary phonon frequency 
luq. This approximation is current-conserving for identi- 
cal molecule- lead couplings = Tr = T and supercon- 
ductor gaps Al = Aft = A. 



Our main results are as follows, (i) For eV S> A, 
the coupling A tends to enhance the excess current I exc , 
which is defined as the difference in current for V — > oo 
when changing normal into superconducting leads, (ii) 
In the subgap regime eV < 2A, the phonon gives rise 
to a quite complicated structure in the I-V characteris- 
tics around each MAR onset at eV = 2A/n (integer n), 
with a pronounced even-odd parity dependence. These 
results can be qualitatively understood within a MAR 
ladder picture. Such a picture has previously been used 
for junctions with a resonant level |2ll| and is here ex- 
tended to include phonon-induced transitions (inelastic 
MAR). Rich features in the I-V curve appear already for 
eV -C hujo, in contrast to normal leads where phonon sig- 
natures (e.g. sidebands) emerge only at eV > TujJq 0,0- 
(iii) For V = 0, we give analytical results for the Joseph- 
son current in the adiabatic limit, fuv <C A <C I\ We 
find a reduction (but no destruction) of the critical cur- 
rent and a changed current-phase relation. These find- 
ings are in qualitative agreement with Rcf . [22( , where the 
opposite limit r <C A has been studied by lowest-order 
perturbation theory in the lead-molecule hopping. 



The outline of the paper is as follows. In Sec.[n] we dis- 
cuss the model and the Keldysh Green function approach 
taken in this work. In Sec. lIIII we present our results and 
provide a physical interpretation in terms of a MAR lad- 
der picture. How such a scattering-type approach can be 
connected to the Keldysh approach is explained in detail 
in the Appendix. Finally, some conclusions are offered in 
Sec. IIVI We put e = h = 1 in intermediate steps. 




II. MODEL AND KELDYSH APPROACH 
A. Model 

We now wish to formulate and analyze a tractable 
model describing the relevant physics of a molecule sand- 
wiched between superconducting leads. Writing the 
Hamiltonian as 

H = uotfb + ]T(eo + XX)4d a + H L + H R + H T , (1) 

we consider one relevant molecular level associated with 
the fermion operator d a for spin a =f , J. and located at 
the energy eq. In Eq. Q we take a linear coupling be- 
tween the molecular charge and the phonon displacement 
X = b + &t , where the boson operator b annihilates a 
phonon excitation. For a justification of this form and 
possible other couplings, see Ref. Q. The leads are de- 
scribed by a pair of standard s-wave BCS Hamiltonians. 
Using the Nambu vector \&J fc = (V'j.fc.T i i) f° r e ^ ec ~ 
trons in lead j — L/ R, we have 

with single-particle energy = k 2 /2m — cf; the 2x2 
matrix acts in Nambu space. In the following, standard 
Pauli matrices in Nambu space are used and denoted 
by u x ,y.z- Using the Nambu vector d = (c£|,g?J) t and 
T = 7rz/ |<o| 2 for (normal) lead density of states Uq, the 
lead-molecule coupling is 

H T = t & jk <rze ±u * vt/2 d+h.c, (3) 

k,j=L/R=± 

where the voltage V enters via the time-dependent phase. 
As we are mostly interested in the V / case, for sim- 
plicity, we consider A > to be real- valued. (In the 
study of the Josephson current, the residual V = phase 
difference across the molecule is of course taken into ac- 
count.) 

Note that Eqs. (J2J and Q assume symmetric coupling 
and identical superconducting gaps. Our approximation 
scheme below will yield a current-conservating result only 
then. In fact, since the calculation of MAR-dominated 
transport is already involved for A = 0, a nontrivial 
current-conserving self-consistent approach covering the 
large transmission limit seems out of question. Note 
that self-consistency is (usually) sufficient to ensure cur- 
rent conservation |23j. Below we instead proceed using 
a perturbation theory with the small expansion parame- 
ter A/T. Under such an approach, current conservation 
is known to only hold for the electron-hole symmetric 
case [24). i.e., eo = 0, with symmetrically arranged leads 
(A L = A R = A,T L = T_r = r). This case is taken 
in what follows. We note in passing that this impor- 
tant issue (and also the tadpole diagram in Fig. was 
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FIG. 1: Self energy due to the presence of the phonon mode: 
(a) 'Fock' and (b) 'tadpole' diagram. The polarization bubble 
(c) leads to the dressed phonon propagator T> (wiggly lines). 
Arrowed lines denote Qq. 

overlooked in Ref. , where the same approximation as 
ours was implemented for normal leads but also used for 
asymmetric cases. 



B. Keldysh approach 

To compute the current- voltage characteristics, we now 
employ the Keldysh Green function technique. This 
method has become a standard approach by now; for a 
review, see, e.g., Ref. [2j|. The Keldysh Green function 
for the d fermion is defined as 

g^(t,t') = -i(f c [d a (t s )dl(t' s ,)}), (4) 

where a, a' (s, s' = ±) are Nambu (Keldysh) indices, 
and Tq is the time-ordering operator along the familiar 
Keldysh contour C [23. Accordingly, t s denotes a time 
taken on branch s of the Keldysh contour. Similarly, we 
define a phonon Keldysh Green function V ss (t, if) for the 
quantity X = b + b^ coupling to the fermion d in Eq. iff]). 
In what follows, we use the check notation (Q) in order 
to schematically indicate the Keldysh structure. In the 
case of electron Green functions, this also includes the 
Nambu structure. Denoting the respective functions for 
A = by Go and 2?o, and using the self-energy diagrams 
in Fig. ^ the dressed Green functions used in our per- 
turbativc approximation follow from the Dyson equation. 
As we consider only the case eo = 0, it can be checked 
that the 'tadpole' diagram does not contribute. 

For the superconducting problem of interest here, it is 
convenient to use the double Fourier representation 

+00 . , 

^')= E / ^-- t ^ t 'e„ m ( w ), ( 5 ) 

n,m= — 00 ^ 

and likewise for all other Green functions and self ener- 
gies. Here we use 

uj n = lu + nV (6) 
for u> within the 'fundamental' domain F defined as 

F = [-V/2,V/2]. (7) 
For fixed w G F , the Dyson equations then lead to the 
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matrix equations 

^nmM = S nm U n f z -T ^ 7j>m(w) (8) 

3=L/R 

Q~ 1 =Qo 1 -% T>- 1 =Vq 1 -ft, (9) 
where the Pauli matrix f 2 acts in Kcldysh space, and 



2^0 



(10) 



The self energy due to tracing out the respective lead is 
given by the Nambu matrix 

^=Va=±,nmM = (11) 

" II HI 

X{oj n t V/2) 5 YK T V/2) 
<Wi ? K ± V/2) <W *(w„ ± U/2) 

with Keldysh matrices Y(w) = -AI(w)/w and 

if z , M < A 



*H = <( iH f 2/ u - 1 -2/ w 
V^A? { 2/_ w 2/„ - 1 



, \w\ > A 



where / w = 1/(1 + e w / feBT ) is the Fermi function. 

Figure n yields for the polarization n and the self en- 
ergy E the following explicit expressions: 

f&>) = -u 2 tr]T | ^ {f K g G , n , m ,{J)f K ) ss ' 

n'm' ^ 

X C»'-™,n'-nK- W )' ( 12 ) 

= *a 2 ^| ^p-;„,„, (w -c') 

n'm' ^ 

X (Tif^0;n'm'(w') f Ar) SS . (I 3 ) 

where 'tr' extends over Nambu space only and tk = f z a z . 
Should the difference u> — to' appearing in Eqs. (|12f> and 
(|13[1 fall outside the fundamental domain F, one has to 
fold it back to F. This is implicitly understood in the 
above equations. 

The steady-state d.c. current through the left/right 
junction then follows in the form 

2~tr [<7zlL/R,nm{u)Gmn(w)\ 

- 

(14) 

This relation constitutes a generalization of the Meir- 
Wingreen formula (2(| to the case of superconducting 
leads. We also define the phonon contribution 



SI ph = 7(A) - J(A = 0). 



(15) 



Current conservation, II = Ir = I, can be explicitly 
verified as follows. In the particle-hole symmetric case, 
current conservation requires that 



This relation is indeed obeyed since the approximate E 
in Eq. (|T3J) also fulfills Eq. JTSJ. 

Using Eq. (|14fl . we evaluate the I-V characteristics for 
A = 0.15r (unless noted otherwise) and fe^T/A = 0.01. 
We truncate the summations such that \uj n \ < uj c — 20A; 
further increase of the bandwidth io c did not change re- 
sults. In practice, the domain F in Eq. 1(7} must be dis- 
cretized. Typically, we found 5u> — 0.008A to be suffi- 
cient for convergence. The matrix inversions in Eqs. © 
and © are then done for each w£f separately, involv- 
ing matrix dimensions of the order uj c /eV . For very small 
eV/A, this becomes quite costly, and we limit ourselves 
to eV/A > 0.15 in the following. 

Fortunately, there are several nontrivial tests that we 
can use to check the scheme. For A = 0, our a ppr oach 
quantitatively reproduces the results of Refs. [l7Lll9Ll2ll |. 
For A = 0, we recover results of Ref. Q when applicable 
(i.e., for e n = 0). As additional check, Green function 
sum rules [25j were verified, such as 



tx{T K g (t,t)) =o 



(17) 



(16) 



at coinciding times. While such relations must hold for 
the exact Green function, it is reassuring to verify that 
our approximation does not lead to violations of sum 
rules. 

Besides the current, we have also monitored the av- 
erage phonon number N p h = (b^b) and the frequency- 
dependent phonon distribution function. We find N p h < 
1 for A/r = 0.15, in accordance with our assumption of 
weak electron-phonon coupling. The phonon distribution 
function revealed renormalizations of the peak position 
away from ujq by a few percent. 

To conclude this section, we note that it is possible to 
establish a close connection between this Keldysh Green 
function approach and a Landauer scattering approach 
incorporating inelastic transitions. Such an 'inelastic 
MAR' picture will in fact be essential in interpreting our 
numerical results in Sec. IIII Al This connection is de- 
tailed in the Appendix. 



III. RESULTS 
A. Subgap regime: Inelastic MAR 

Let us then turn to results for the I-V curve. We 
start with the subgap regime, where MAR provides 
the dominant transport mechanism. In particular, for 
2A/(n + 1) < eV < 2A/n (n integer), there is a total 
number n of Andrccv reflections for electrons or holes 
within the superconducting gap. The I-V curve for 
r = 2A and uiq = 0.2A is given in Fig. [21 where SI p h 
is always negative. Note that in this fully transmitting 
limit, the I-V curve for A = is smooth and does not ex- 
hibit the MAR 'cusps' encountered at lower transmission 
[T?| . Phonons now restore such features near MAR on- 
sets, with pronounced even-odd 'parity' effects: For even 
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FIG. 2: Phonon difference current 1151 for hu)o = 0.2A and 
r = 2A. In all figures, currents are given in units of eA/ (2nh), 
and dotted lines are guides to the eye only. Left inset: Same 
as function of 2 A/el/. Right inset: Part of the total I-V 
curve (note the scales) . 




FIG. 3: MAR ladder picture with phonon emission. Here 
eV is slightly below A: for an electron incoming from the left 
side, we have one hole (open circle) and two electron (filled 
circle) segments. Dashed lines indicate possible trajectories 
after single phonon emission involving either hole or electron 
segments. There is also a MAR path (not shown) for a hole 
entering from the right side, with one electron and two hole 
segments. 




FIG. 4: Same as Fig.|5|but for huio = 1.8A. The upper inset 
gives SIph as a function of 2A/eV, the lower inset gives the 
low- voltage part of the total current. 



(odd) n, SIph shows valleys (peaks) around eV = 2A/n. 
This is clearly seen in the left inset of Fig. [21 for n up 
to 12 [23. Note that the dip at n — 8 is somewhat 
shifted, presumably due to phonon renormalization ef- 
fects. However, the appearance of even-odd parity oscil- 
lations is quite distinct and surprisingly regular given the 
complexity of this system. 

In order to achieve a physical understanding of this 
even-odd effect, it is useful to invoke a MAR ladder 
picture in energy space, including inelastic transitions 
caused by phonon emission. A schematic description of 
the MAR ladder picture in this limit is given in Fig. [3] 
Superconductor spectra are placed next to each other, 
but because of the presence of the bias, electrons emit- 
ted from left to right gain an energy eV in this process. 
Holes that are reflected, traveling from right to left, also 
gain this energy. 

Usuallvthe scattering approach is used to develop this 
picture [2lj . but its straightforward implementation en- 
counters conceptual difficulties in the presence of inelas- 
tic phonon transitions. Fortunately, a different route re- 
solving these difficulties can be formulated via the above 
Keldysh approach, whose formal justification is detailed 
in the Appendix. Here electron and hole propagators in 
energy space are coupled to each other through suitable 
Andreev reflection matching conditions and the phonon 
self energy S. For a qualitative explanation of the even- 
odd effect found in the full calculation, cp. Fig. it is 
sufficient to restrict the full MAR ladder picture to single 
phonon emission processes, where electrons (holes) lose 
(gain) the energy Huq. Since N p h < 1, emission dom- 
inates over absorption and multi-phonon processes are 
rare. 

In Fig. [3J the two superconductors are positioned at 
the same chemical potential, but electrons (from left to 
right) and holes (from right to left) 'climb' the MAR lad- 
der by gaining eV for each Andreev reflection. The higher 
the total number of Andreev reflections in one cycle, the 
larger the total charge transmitted. Since we consider the 
high transmission limit where high-order MAR processes 
are not penalized, the current is therefore expected to 
increase (decrease) if phonon emission is able to increase 
(decrease) the number of Andreev reflections in a MAR 
cycle. For eV slightly below 2A/n with even n, we then 
argue as follows; for n = 2, see Fig. [21 For small en- 
ergy transfer Hu>o, if a phonon is emitted during an elec- 
tron segment, MAR trajectories in energy space are not 
drastically modified in the sense that the number of An- 
dreev reflections stays unaffected. However, if a phonon 
transition occurs during a hole segment, the MAR lad- 
der is shifted upwards by Hujo and the last hole on the 
MAR ladder can be scattered into the continuum (left 
electrode in Fig. [3J instead of being Andreev reflected. 
Consequently, one Andreev reflection is lost and hence 
the current is expected to decrease. This argument ap- 
plies both to incoming electrons and holes, and explains 
why current valleys are observed for eV ~ 2A/n with 
even n in Fig. [21 On the other hand, consider eV slightly 
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B. Excess and Josephson current 




FIG. 5: Same as Fig.^but for F = 0.5 A. The left inset gives 
SIph versus 2A/eV, the right inset the total current. 



above 2A/n with odd n. Reiterating the above analysis, 
now phonon emission during a hole segment tends not to 
affect the number of Andrccv reflections. If the phonon 
is emitted during an electron segment, however, an ad- 
ditional Andreev reflection has to take place to complete 
the MAR cycle, leading to a current peak for eV m 2A/n 
with odd n. 

Next consider Hcoo = 1.8A but otherwise identical pa- 
rameters, see Fig. 0J where 5I p h can be positive and again 
shows oscillations near the MAR onsets, which are less 
pronounced for small n = 2A/eV. Remarkably, even for 
small voltages, eV <C hwo, a rather complicated subgap 
structure is caused by the phonon. At such low voltages, 
this is only possible via MAR, for otherwise electrons 
or holes do not have enough energy to emit a phonon. 
The broad minimum corresponding to n = 2 observed 
in Fig. |2 has now vanished. Let us invoke the MAR 
ladder picture to rationalize this effect. For eV < A, 
by emitting a high-energy phonon (Hluq > A), the last 
electron on the MAR ladder can now be scattered back 
inside the gap instead of heading to the continuum. This 
increases the number of reflections and thus the current. 
As a phonon emitted during the hole segment has the op- 
posite effect (see above), the net outcome of the higher 
phonon frequency is to suppress the valley at eV ~ A 
expected for small ljq. Figure 0] also shows a dip in the 
current at eV « 1.8 A, representing a phonon backscat- 
tcring feature at eV = Hloq. The scaling of this dip with 
u)q was confirmed by additional calculations. One can 
also see a two- phonon feature at eV = ha>o/2 in Fig. 0] 

Let us then briefly go back to Hujq = 0.2 A (cp. Fig. [2J . 
but now for V = 0.5A, see Fig. [3J For small T, quasi- 
resonances appear plj and cause additional features, 
e.g. negative differential conductance portions in the I- 
V curve. The MAR ladder picture then has to include 
both quasi- resonances (cp. Ref. pl|) and phonon tran- 
sitions, which is possible but beyond the scope of this 
paper. Fig. [5] shows that the even-odd parity effect re- 
quires large T to be observable. 



Finally, we briefly discuss the limits of very large 
and/or zero voltage. Starting with the first case (ex- 
cess current), we have computed the difference 5I e xc,ph 
between the excess currents l eX c with and without the 
phonon. This calculation has been done at eV = 10A, 
which according to the discussion in Ref. |2l| is certainly 
large enough for our parameters. For the case of high 
transmission encountered here, we find that phonons gen- 
erally enhance the excess current. To give a concrete 
example, for Hojo = 0.8A, V = 2A and A = 0.5A, we 
find SI exCjP h/I e xc ~ 0.07. A similar current enhancement 
at high transmission was also found for environmental 
Coulomb blockade in superconducting jun ctions, and has 
been explained as 'antiblockadc' effect [28j. As such, this 
effect of the phonon may not be too surprising. 

Second, the equilibrium Josephson current has been 
calculated by adopting our approach to the Matsubara 
representation. For arbitrary parameters, it is straight- 
forward to numerically compute the full current-phase 
relation. However, the current-phase relation in the adi- 
abatic phonon regime defined by 

huj « A « r (18) 

can even be calculated analytically, with the result 

7(0) = (eA 2 /2h) T sin(0)/£ o (0), (19) 

where 

E (<f>) = A[l - T sin 2 W-/2)] 1 / 2 (20) 

is an Andreev bound state energy [Tt| in the junction 
with an effective transparency 



To 



i + A 2 /4r 2 



(21) 



The ^-dependent shift without any broadening of the An- 
dreev level caused by the coupling to a phonon mode is 
characteristic for the coherent limit and decreases the 
critical current. Very similar results were reported in 
Ref. in the opposite limit T < A. 



IV. CONCLUSIONS 

In this work, we have for the first time theoreti- 
cally explored noncquilibrium molecular transport with 
superconducting leads in the coherent phonon regime. 
Phonons reveal a rich subgap structure even for voltages 
well below the phonon frequency, including a surprising 
even-odd parity effect near the MAR onsets. This effect 
can be largely understood within a physically appeal- 
ing inelastic MAR ladder picture, based on the assump- 
tion that single phonon emission processes dominate. For 
stronger clcctron-phonon couplings, also multiple phonon 
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processes and/or absorption becomes important, and the 
practical usefulness of such a scattering approach is less 
clear, although it can be formally derived as outlined in 
the Appendix. 

Let us also offer a brief outlook. Besides the obvi- 
ous interest to experimentally probe the effects predicted 
here, the problems raised above deserve further theoret- 
ical work. In the case of strong electron-phonon cou- 
pling, different approximations will be necessary. More- 
over, other quantities such as shot noise or frequency- 
dependent noise [2j| deserve attention. Noise can yield 
information about the effective charge involved in the 
transfer process, and how this charge is affected by 
phonon transitions remains to be explored. It would also 
be interesting to formalize the notion of inelastic MAR 
spectroscopy, potentially allowing to infer the electronic 
structure of molecules from the superconducting current- 
voltage characteristics. Such approaches have already 
been very useful in nanoscale break junctions, and are of 
obvious interest in the molecular electronics context as 
well. It is safe to conclude that superconducting molecu- 
lar transport definitely warrants further surprises in the 
future. 
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(r) (r) 

where = . Here we consider the simplest 



J, — J— I* 3 k ----- 

possible situation, where phonon rcnormalization effects 
like the polarization bubble in Fig. ^ are completely dis- 
regarded. In that case it is sufficient to simply take an 
equilibrium average over the phonon subsystem, while 
otherwise one should also take into account the equation 
of motion for X(t). 

The solution of Eq. adapted to ^j(t), can 

be written as the sum of a free (homogeneous) solu- 
tion tyjQ (t) , describing an incoming quasiparticle in lead 
j = L/R, plus a scattered part due to the interaction 
with the molecular level, 



* r, (t) = * *o dt> gf\t-t>) (24) 



xa z e 



±i^Vt'/2 d ^ 



which is expressed in terms of the retarded Green func- 

(r) 

tion g\ for uncoupled electrodes whose Fourier trans- 
form is given by 



9l KM) - 9r \u) - -j- 



A(7 T 



V(w + iOp - A : 



When summed over r, this essentially yields the retarded 
components of jj in Sec. [H] 

In what follows, we use s = {1,2,3,4} in order to 
label scattering processes corresponding to electron- or 
hole-like quasiparticles incoming from the left (s = 1,2) 
or right (s = 3, 4) electrode. In Fourier representa- 
tion, where E is the energy of the incoming quasiparticle 
(\E\ > A) and E n = E + nV, we write 



INELASTIC MAR: SCATTERING APPROACH 

In this appendix, we outline the derivation of the Lan- 
dauer scattering approach to inelastic MAR incorporat- 
ing phonon transitions. Such an approach is formulated 
here using the equation-of-motion method, where trans- 
fer matrices are expressed in terms of Keldysh Green 
functions. This provides a formal justification for the 
intuitive MAR ladder picture used in Sec. IIII Al 

Using the Hamiltonian we can write down the 
equations of motion for the quasiclassical envelope func- 
tion describing right/left-moving quasiparticles (r = 
±) in the left or right lead (j = L/R = ±) with momen- 
tum rkp + k. Here |fc| <C kp is assumed, and with the 
Fermi velocity v F , we obtain 

(id t - rkv F a z - Aoi)*5=± fc (t) = t a z e ±l ^ vt ' 2 d(t). 

(22) 

Similarly, the equation of motion for the molecular 
fermion (Nambu spinor d) is 

(id t -XX(t)a z )d(t)=t a z J2 e^ vt/2 ^\t), (23) 

j=±,r 




(25) 



and introduce electron- and hole-type spinors, 



<!>'• 



h 



h — 



These spinors are defined on a chiral space, with the two 
entries corresponding to the right- and left- moving parts. 
From Eq. (|24p. we proceed to derive a first set of Andreev 
reflection matching equations for , 

<f> e ( a ™ \<b h -LX I-,, 2 „ 2 \ ( S ^/ u e \ 



a n J Rn + n0 ^ E ~ E > { 8 s3 /u E 



where a n = (E n — yf (E„ + iO) 2 — A 2 )/ A. Furthermore, 
the Ue,Ve denote the n = entries in Eq. Note 
that these matching equations are not modified by the 
presence of phonons, see Ref. [30l |. 
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Iterating Eq. I|23[) and averaging over the phonons, we 
find 

(idt - E) o d = t a z e^' vt ' 2 0>f (i), (26) 

j=±,r 

where o is a shorthand for convolution. The phonon self 
energy entering Eq. I|26[) is given by 



E(M') = i\ 2 V>(t,t')<j z g«(t,t')a 



(27) 



Note that this ignores the polarization-bubble renormal- 
ization, see above. For the four different scattering pro- 
cesses indexed by s, it is convenient to introduce matrices 
S nm according to 



E 



(JlC 

2^ 



—iuj n t-\-iuj m t 



E nm (w) (28) 



with a s = 1 (a s = -1) for s e {1,2} (s e {3,4}). To 
keep the notation simple, the s-dependence of S nm is not 
exhibited explicitly. 

Using Eq. H26J1 . wc then obtain a second set of matching 
equations. For s = {1,2}, 

L?,n+1 \ M nJ Ln ' / , V nm ^ Lm ' ^nm ^ Lm } 

n l 

$n 1 = (T h V 1 ^'r + "V (S^ $ r "h^ 6 $f ) 

^R,n— 1 V n/ ^in ' / y V mn Lm ^ '-'nm ^ Lm) i 

m 

while for s = {3,4}, 

L,n— 1 n— 1 -Rn 1 / . \ nm Km ' ^nm ^Rm) 



®L,n+l — T.n +1 <& Hn + ^ (<5ram ^flm + ^nm ^ 



with 



= -1/(1 T *^ / V2r), r^V^ A - ±^ /h /2r, 



where EfJ h = E + (n ± 1/2) V". Furthermore, transfer 
matrices linked to phonon-induccd transitions are given 
in terms of the phonon self-energy matrix elements intro- 
duced in Eq. (p^l 



gee geh 
g/ie chh 



(30) 



-£ n (£) £ 12 (E)\ 



where the Pauli matrices Tj, jZ operate in the space of the 
Q e ' h spinors and the a indices in T, aa refer to Nambu 
space. Note again that these quantities all carry an im- 
plicit s-dependence. 

Finally, after straightforward but somewhat tedious al- 
gebra, the d.c. current is expressed in terms of the above 
scattering states, 



/ = evp 



dE- 



\E\ 



|B|>A v^ 2 - A 2 



/*EE (3D 



E [(^+^)<(£)r^| n (£) 



a—e/h 



+ (S sS + d s ^li(E)T^ a Ln (E) 



Rmj 



Here, electron/hole transfer matrices are given by 



T e /h= ( l/t e J h (r*J h /fJ h 
V r e n /h /t e n /h l/f n /h * 



Similarly, the full Keldysh Green function Q can also be 
constructed from the complete set of scattering states 



(29) <^(E). 
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